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ABSTRACT 

Observations show that propagating magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmo- 
sphere. The technique of MHD seismology uses the wave observations combined with MHD wave theory to 
indirectly infer physical parameters of the solar atmospheric plasma and magnetic field. Here we present an an- 
alytical seismological inversion scheme for propagating MHD waves. This scheme uses in a consistent manner 
the observational information on wavelengths and damping lengths, along with observed values of periods or 
phase velocities, and is based on approximate asymptotic expressions for the theoretical values of wavelengths 
and damping lengths. The applicability of the inversion scheme is discussed and an example is given. 
Subject headings: Sun: oscillations — Sun: corona — Sun: magnetic fields — waves 



1. INTRODUCTION 

Standing transverse magnetohydrodynamic (MHD) 
waves in coronal loops were first detected in 1999 by 
lAschwanden et all (1999) and iNakariakov et all (119991) in 
observations made with the EUV telescope on board of the 
Transition Region and Coronal Explorer (TRACE). Since 
then the detection of these standing MHD waves has been 
confirmed and in addition damped standing MHD waves 
have been observed in hot coronal loops by the SUMER 
instrument on board SOHO More recent observations of 
transverse loop oscillations with high spatial and temporal 
resolutions have been made with instruments on board of 
STEREO and SDO (see e.g.. IVerwichte et al.1 120091 1 201 0[ 
Aschwanden & Schriiven 1201 1[ IWhite & Verwichtef l2012t 
Wang et al.l 120121) " The transverse loop oscillations have 
periods, T, of the order of 2-10 minutes and compar- 
atively short damping times, tq, of the order of ^ 3-20 
minutes. There is general consensus that the transverse loop 
oscillations are standing kink MHD waves. Kink refers to 
the azimuthal wavenumber m — 1 in a system of cylindrical 
coordinates with its z-axis along the axis of the loop. The 
MHD waves have to be kink because only for this value of 
the azimuthal wavenumber the axis of the loop and the loop 
as a whole are displaced. 

A possible explana tion of the observed rapid d amping is 
resonant absorption. iRuderman & Roberts! (12002b were the 
first to suggest that the observed rapid damping of the trans- 
verse oscillations of coronal loops could be explained by reso- 
nant absorption. In the c ontext of the heating of solar plasmas 
iHollweg & Yang! (fl988) have predicted that oscillations in 
cor onal loops are to under go rapid damping. In the same con- 
text |GooiiemeraL] d 19921) derived analytical expressions for 
the frequency and the damping rate of quasi-modes in static 
and st ationary equilibrium models. Ruderma n & Roberts! 
j2002h focused on proving the principle of resonant absorp- 
tion as damping mechanism for the transverse standing MHD 
waves in coronal loops and con sidered one specific numerical 
example. Goos sens et al.l d2002l) looked at the damping times 
of 1 1 loop oscillation events and basically confirmed that res- 
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onant abso rption can explain the obse r ved damping a s sug- 
ges ted by IRuderman & Robertsl (120021) . I Arregui et all (2007) 
and lGoossens et all (|2008 ) showed that resonant absorption of 
kink MHD waves can explain the observed periods and damp- 
in g times. 

iGoossens et al.l (120091) showed that long-wavelength kink 
MHD waves are primarily driven by magnetic tension force 
and hence their behavior is more Alfvenic than fast mag- 
netosonic. The condition of long wavelegth is that the 
wavelength is much longer than the radius of the tube. 
IGoossens et al.l d2012l) showed that kink MHD waves prop- 
agate vorticity and that the fundamental radial modes of kink 
waves are surface Alfven waves. Note t hat the term A lfvenic 
was alre a dy int roduced in the 1970s by Ionson ( 1978) and by 
IWentzell (11979). The surface Alfven wa ve chara c ter of kink 
MHD waves was already pointed out bv lWentzefl d 19791) . An 
important property of surface Alfven waves is that in a non- 
uniform plasma their frequency is in the Alfven continuum so 
that they underg o continuum damping or resonant abso rption 
(see, e.g., Hase gawa & Uberoil fl982: Go edbloedl 19831) . 

Standing MHD waves are relatively rare phenomena since 
they need an energetic event like a solar flare or strong vor- 
tex shedding in order to get excited. A more recent devel- 
opment over the last five years is that observations show that 
small amplitude propagating transverse MHD waves are al- 
most everywhere in the solar atmosphere. They were first ob- 
served with the Coronal Multi - chann el Polarimeter (CoMP) 
instrument by ITomczvk et al.l d2007l) and subsequently by 
iTomczvk & Mclntoshl (120091) and others. The physical mech- 
anism of resonant absorption does not make a distinction be- 
tween standing MHD waves and propagating MHD waves. 
The physical mechanism is the same and the mathematical 
analysis is largely the same. For a standing MHD wave the 
axial wavenumber k, is specified and the frequency a> is deter- 
mined. For a propagating MHD wave the frequency a> is spec- 
ified and the wavenumber k z is computed. The damping time 
td for standing waves is replaced with the damping length 
Lo- In what follows we shall present an inversion scheme 
for propagating MHD waves that uses observed values of the 
wavelength and damping length in combination with asymp- 
totic analytic expressions to compute theoretical values of the 
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wavelength and the damping length. The theory of resonant 
damping for p r opagati ng MHD waves has been developed by 
Terradas et al. (2010b) and used to explain CoMP observa- 
tions by lVerth et al.l (120101). Subseque nt extensions to the the- 
ory were made by lSoler et al.l (1201 lallbl lcl). 

The present investigation is inspired by earlier work on 
seismic inver s ion for standing MHD waves in coronal loops. 
lArregui et al.l (120071) were the first to carry out a seismo- 
logical investigation of standing transverse MHD waves of 
coronal loops that used the observed values of both the pe- 
riods and the damping times in a consistent way. Their seis- 
mic inversion was fully numerical and the theoretical values 
of the periods and damping times were obtained by numer- 
ical eigenvalue computations for one dimensional cylindri- 
cal equilibrium models. As a result the inversion scheme 
could turn out to be rather invo lved. In a subsequent in- 
vestigation iGoossens et al.l (120081) used approximate asymp- 
totic expressions for the period and damping rate to work out 
an analytical seismological inversion scheme. This analyti- 
cal inversion scheme has the big advantage that it is amaz- 
in gly simple to u s e whe n compared to the numerical scheme 
of lArregui et al.l (120071) . The downside of this scheme is 
that its simplicity might invite us to use it outside the do- 
main of validity of the approximate expressions for the pe- 
riod and damping rate. Analytical schemes are also very 
useful for parameter inference in the Bayesian framework. 
The main potential of Bayesian inference is the consistent so- 
lution to the inverse problem using the forward model and 
the available observational information. The unknown pa- 
rameters are constrained by data and uncertainties are cor- 
rectly propagated from observed data to inferred parame- 
ters. Analytical forward and inverse problems great l y sim- 
plify this task, a shown by lArregui & Asensio Ramosl d201 ll) 
using the analytical scheme of IGoossens et al.l (120081) to con- 
strain coronal loop properties. In addition, the inversion 
scheme has been used beyond the context of coronal loop 
oscillations. For example, it has been used to perform 
seismology of thin threads of solar prominences, since the 
transverse oscillations of both coronal loops and prominence 
threads can be studied using the same physical model (see 
details in, e.g . ISoler et all 120101: lArregui & Bal lesteri |20TT1: 
lArregui et al.ll2012l) ISoler et all (1201 lbl) suggested that the 
inversion scheme of IGoossens et alj (120081) could be adapted 
to the case of propagating waves. Here we fully explore the 
seismic inversion for propagating MHD waves. 

2. ASYMPTOTIC ANALYTIC EXPRESSIONS FOR THE 
WAVELENGTH. 

In this section we derive an expression for the wavelength 
of a propagating wave with a given frequency in the thin tube 
approximation. In the following section we shall derive an 
expression for the damping length in the thin tube approxi- 
mation and thin boundary approximati on. Part of the infor - 
mation in this section is presented in lTerradas et alJ (120 1 Obi) . 
iTerradas et al.l (1201 Obi) investigated the damping by resonant 
absorption of propagating MHD waves. They obtained simple 
expressions for the wavelength and the damping length in the 
thin tube approximation and the thin boundary approximation. 
In addition they went beyond the thin tube and thin boundary 
approximation by using numerical resistive calculations. 

The analytical expression that we shall use for the wave- 
length is obtained by (i) adopting the thin tube (TT) approxi- 
mation for MHD waves and by (ii) modeling the wave guide 
as a uniform cylinder with a straight magnetic field along 



the z-axis. For a standing MHD wave the axial wavelength 
is specified by the dimension of the flux tube and the cor- 
responding frequency is determined by the dispersion rela- 
tion. In that case the thin tube approximation means that axial 
wavelength is much longer than the radius R of the tube so 
that k z R « 1. In case of propagating MHD waves the fre- 
quency is specified and the wavelength is determined by the 
dispersion relation. The thin tube approximation now means 
that during one period specified by the frequency a> a signal 
traveling at the Alfven speed Va can cross the wave guide in 
the radial direction many times or 



v A /R 



« 1. 



(1) 



The TT approximation means the wavelength is independent 
of the radius and that effects due to non-zero radius are absent 
as far as the wavelength is concerned. The choice of a uniform 
equilibrium model means that effects due to stratification are 
absent. The wave guide is modeled as a cylindrical plasma 
with constant density p; embedded in an external plasma with 
constant density p e . The wave guide is basically a density 
enhancement with p e < Pi- The magnetic field is constant and 
has the same strength both inside and outside the tube. 

Our starting point is the well-known thin tube approxima- 
tion of the dispersion relation for non-axisymmetric MHD 
waves on a uniform cylinder with a straight magnetic field 
along the z-axis, namely 



Pi(w 2 - W A i ) +p e (0) 2 



W i,e) 



0. 



(2) 



The subscripts 'i' and 'e' refer to quantities respectively in the 
tube and in the external plasma surrounding the tube, p is the 
density, wa is the local Alfven frequency and va is the local 
Alfven velocity. They are defined as 



w A = k z v A , 



B 



VA 



(3) 



where k z is the longitudinal wavenumber, B is the magnetic 
field strength and /j.o is the magnetic permeability of free 
space. Note that the dispersion relation (Equation (f2]i) is inde- 
pendent of the azimuthal wavenumber m so that Equation (O 
applies to all non-axisymmetric MHD waves. Our interest is 
focused on kink waves which have m = 1 . 

For standing waves k z is specified and the dispersion rela- 
tion Equation © is solved for the frequency a>. For propa- 
gating waves we consider waves that are generated at a given 
position with a real frequency and we have to solve Equa- 
tion (f2]l for k z . The result is 



, Pi +Pe _ , 



(4) 



Let us convert frequencies u to periods T = 27r/a> and 
wavenumbers k z to wavelengths A = 2n/k z and rewrite Equa- 
tion © as 

V2v Ai r 

A = ' . (5) 

MO 

In Equation (|5]l £ is the density contrast, namely 



(6) 



The function A (£) is defined in Equation (4) of lGoossens et al.l 
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d2008h ; 



A(0 = 



ill 



1/2 



(7) 



Equation (0 is our first key equation. It expresses a relation 
of the wavelength A, and the period, T, which are two observ- 
able quantities, in terms of the internal Alfven velocity, VA,i 
and the density contrast, (, which are two quantities that we 
aim to determine by seismic inversion. Let us recall that Equa- 
tion (0 has been obtained by use of the TT approximation for 
a uniform cylindrical tube. Effects from non-zero radius and 
stratification are absent from Equation (0. 

Let us now look at Equation (0 from the seismic point of 
view. If we have observed values of the wavelength, A, and 
the period, T, and we convince ourselves that Equation <(5j is 
a good first analytical approximation of the wavelength then 
we can invert Equation (0 for either vaj or Z\ Actually, we 
shall do both. Let us first solve Equation (0 for VA,i. We 
prefer to use dimensionless quantities and so we introduce y+ 



y* = 



VAjT 



From Equation (0 we obtain 



y* 



V2 



A(0- 



(8) 



(9) 



Since we have not any information on f, it might appear that 
Equation (0 is not very helpful. However, closer inspection 
reveals that it contains valuable information. First note that 
for given observed wavelength A, Equation (0 is a parametric 
representation of (or equivalently of v>A,i) m terms of (. In 
order to stress this point we define the function F\ iis by use of 
the right hand member of Equation (0 as 



V2 



(10) 



Second, note that the function F i_* is strictly decreasing and 
in addition note that 



max(Fi 1 +)=Fi f +(l) = l, lim Fi,*(£) = — . (11) 



This means that 



-p< y+ < l, 

1 A A 

< Va i < — . 

V2 T a " T 



(12) 



(13) 



Inequality (fT~3T > tells us that the Alfven velocity is constrained 
to a narrow range. Whatever the density contrast is, the 
Alfven velocity is in between A/ ^JlT and A/T. 

Let us stress again that the seismic variable y+ contains in- 
formation of two observables, namely the wavelength, A, and 
the period, T. However, we note that the ratio A/T can be 
expressed as 

| = v ph , (14) 

where v p h is the phase velocity of the propagating wave. Us- 
ing Equation (fl4l the definition of becomes 



y* = 



VA,i 
Vph 



(15) 



Therefore, for practical purposes we only need information 
about one observed quantity, i.e., the phase velocity, v p h. 
Hence Inequality ( TT3l can be rewritten as 



1 



ph 



< VA,i < V 



ph- 



(16) 



With the help of Inequality ( TT2l we can refine the definition 
of F\j, and replace Equation ( TTOb with 



Fi.* : [1, «,[-»]—, 1], ( => Fi,*(0 = 
V2 

Let us now solve Equation (0 for £ and find 

1 



1 (S+i 

V5l i 



M 1 - 



(17) 



(18) 



Equation ( TT~8b is the twin of Equation (0. For given observed 
wavelength, A, and period, T, or alternatively for a given ob- 
served phase velocity, v p h, Equation ( fT8l is a parametric rep- 
resentation of f in terms of (or equivalently in terms of 
VA,i). In order to stress this point we define the function Gi * 
by use of the right hand member of Equation ( fTSI ) as 



Gi,* !] ^ 

V2 



Gi,*Cy*) 



1 



V,-i 



The function Gi,* is strictly decreasing and 

Gi,*(l) = 1, lim Gi^(y) = 00. 

y-»l/V2 



(19) 



(20) 



With this information on the function Gi* we can refine its 
definition (Equation (fT9l)). Combined with the definition of 
F\ (Equation ( fT71 )) we obtain the following prescriptions for 
the functions F\^ and Gi,*, namely 



+ 1 



'l,* 



1 1 (t 

[1, «»[-»[—, 1], f => F u (0 = -= 

V2 V2 \ 4 



]— , 1] — > R, => Gi,*(y*) = — = — -. 
V2 2yi - 1 



(21) 
(22) 



Of course, Gi,* is the inverse function of Fi_*, i.e., Gi,+ = 
Fj" ^ and conversely G\\ = F\ jir . 

Let us recapitulate what seismic information we have de- 
duced from the observations. The quantities f and are in 
the following intervals 



f e/ f = [1, oo[, e 7 V 
and are related to one another by 



1 

Vl' 



y* 



F h J(), f = Gi,*(y+). 



(23) 



(24) 



The functions Fj * and Gi* are defined by Equations d2Tb and 
(l22l . When the wavelength, A, and the period, T, are known 
from observations, or alternatively when the phase velocity, 
Vph, is known, then there is infinite number of couples (<T,;y*) 
that reproduce the observations. We can let f vary over the 
interval [1, ; 00 [ and compute for each value of £ the cor- 
responding value of = Fi *(£), or conversely, let vary 
over the interval [ 1 / v2, 1 [ and compute for each value of y+ 
the corresponding value of £ = Gi,+(y). 
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3. ASYMPTOTIC ANALYTIC EXPRESSIONS FOR THE 
DAMPING LENGTH. 

Let us now turn to the damping length. In order for the kink 
MHD waves to be damped by resonant absorption additional 
physics has t o be introduce d in the equilibrium model. As 
explained by iTerradas et al.l d2010bl) the required additional 
physics is non-uniformity of the local Alfven velocity. For 
a constant magnetic fi eld this implies a non-uniform density. 
ITerradas et al.l (1201 Obi) derived an asymptotic expression for 
the damping length. This asymptotic expression is derived 
in the approximation that the non-uniform layer is thin. This 
is the so-called thin boundary approximation. In what fol- 
lows we shall refer to it as the TB approximation. The true 
density discontinuity is replaced by a continuous variation 
of density in the interval ]R - I R + |[. / is the density 
inhomogeneity length scale. The use of the TB approxima- 
tion results in the mathematical simplification that the MHD 
waves can be solutions for uniform plasmas that are connected 
over the dissipative layer by jump conditions. In this way 
we can avoid solving the non-ideal MHD wave equations. 
The jump condition for the ideal A lfven singularity was in- 
troduced on an intuitive manner by Hollwe g & Yangl (1988) 
and put on a firm math ema tical basis by Sakurai et al. ( 1991), 
iGoossens et al.l d 19951) . and lGoossens & Ruderman| (19951) for 
the driven problem and by iTirrv & Goossensl d 19961) for the 
eigenvalue prob lem. Jump conditions were disc ussed and 
used i n, e.g., IGoossens et alj d2006l 120111) and iGoossensI 
(2008]). In case of the TB approximation combined with the 
TT approximation the inclusion of the effect of a non-uniform 
layer in the dispersion relation is relatively simple. Dispersion 
relation (Equation (f2]i) is now modified as 



Pi(a) 



l.i)+Pe(t 



W A.e) 



m/r A 



-p(r A )|A A | Pi(w2 -<' )Pe(w2 -^ X25) 

In Equation d25l l r A denotes the position of the Alfven reso- 
nance. In the TB approxi mation it is natura l to adopt r A = R 
but see the comment by Sole r et ail (120091) . The use of the 
jump conditions is not restricted to thin non-u niform layers as 
can be seen from, e.g.. lTirry & Goossensl d 19961) . The solution 
of Equation d25l l is now a complex wavenumber k, — £r + iki. 
Its imaginary part k\ re flects the spatial damping of the wave. 
ITerradas et al.l (1201 Obi) solved Equation (l25l l in the assump- 
tion that k\ « &r, i.e. that damping is weak. In that case the 
real part £r is given by Equation (|4j, i.e., £r = The result 
of the analysis for ky is 



k\ n m {pi — p e ) 1 
K ~ 8 R pi + p e |dp/dr| 



(26) 



The damping length Ld is defined as Ld = 1/&I- Equa- 
tion d26l i can be rewritten as equation for the damping length 
over wavelength, namely 



Ld 
A 



4 R Pi + p e 
TT 2 m (p; - p e ) 2 



(27) 



Note that the effect of the inhomogeneous layer is contained 
in the value of the spatial derivative of density at the resonant 
position. Note also Equation (f27]i is valid for all m > 1 . In 
what follows we concentrate on kink waves with m = 1 . It 
makes sense to adopt 

dp Pi - Pe 



dr 



= a- 



I 



(28) 



a is a numerical factor that depends on the profile of the vari- 
ation of density in the inhomogeneous layer. For a sinusoidal 
density profile a = n/2, while for a linear profile a = 1. With 
Equation (l28T i th e expression in Equatio n (l27T i is now (see also 
Equation (13) of ITerradas et"aTll2010bl) 



Ld _ F 1 Pi+Pe 
A II R Pi-p e ' 



(29) 



where F = Aa/n 2 is a numerical factor. F has a in it. For a si- 
nusoidal profile F - 2/n and for a linear profile F — A/n 2 . In 
what fo llows we shall adopt the value F = 2/n for a sinusoidal 
profile. Terrada s et al J (120 10b) have gone beyond the TB and 
the TT approximation in an attempt to assess the accuracy 
of the analytical expression of Equation d29t . Their conclu- 
sion was that Equation (|29| | is accurate far beyond its interval 
of applicability. On another note we would like to point out 
that Equation d29l ) is obtained from an asymptotic analysis of 
the boundary value problem. This asymptotic analysis adopts 
an exponentially decaying solution. In principle this expo- 
nentially decaying solution is only valid for z — > °° in the 
same way as the exponentially decaying solution in time is 
strictly only valid fo r t — > oo. However, numerical analysis by 
Pasco e et al.l d2012l) shows that deviations from the exponen- 
tial damping only occur for relatively small density contrasts 
and short distances. 

At this point it is important to note that in the TT and TB 
approximations the analysis for temporal damping of standing 
waves and spatial damping of propagating waves is equiva- 
lent. The expressions for td/L and L^/A are exactly the same. 
A consequence of this fact is that we can repeat Section 3 of 
iGoossenset al. (2008) with the period T and damping time 
To replaced with wavelength A and damping length Ld. We 
rewrite Equation ( f29b in terms of namely 



A 



2£+l 1 
tt£- lT/R' 



(30) 



This is the second key equation of the present investigation. 
Let us now look at Equation d30b from a seismic point of view. 
If we have observed values of the wavelength A and the damp- 
ing length Ld and we convince ourselves that Equation d30l l 
is a good first analytical approximation of the damping length 
then we can invert Equation (l30l > for either f or l/R. Actually 
we shall do both. Let us first solve Equation (f30b for £ and 
find 

( =2 l-4 ■ (31) 

2R A 1 

The reader must keep in mind that we are looking at this equa- 
tion from a seismic point of view. Hence wavelength and 
damping length can be considered as known from observa- 
tions and it makes sense to denote ttLq/A as a constant C. In 
addition it is convenient to abbreviate 1/2R as x. Hence 



I 

2R' 



C = 



TlLr 



Equation (f3TT > can be rewritten as 

Cx+ 1 



Cx-1 



(32) 



(33) 



Here are several observations to be made. First, since £ > it 
follows from Equation (l33l that for a given ratio Lu/A there 
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is a lower bound for the inhomogeneity length scale, namely 

(34) 



_ I A _ 1 



Second, realize that Equation ( l33l is a parametric representa- 
tion of ( in terms of x — 1/2R. In order to make this point very 
explicit, we introduce the function G 2 ,+ defined as 



G 2 ,* :]i 1] -» 
It is easy that show that 



x => G 2 (x) = (35) 
Cx - 1 



-2C 



(Cjc- l) 2 



<0. 



(36) 



Hence ( is a decreasing function of x = Z/2/?. It attains its 
absolute minimum for x = 1/2R = 1. This minimal value is 



£min - G 2j *(l) 



C+ 1 

c- r 



(37) 



Conversely £ attains its maximal value in the limit x — > 1/C, 
namely 

Urn G 2 ,*(jc) = 00. (38) 

With this information on the function G 2 * we can refine its 
definition (Equation (l35l l) as follows 

1 C + 1 Cx + 1 

G 2> * :]-, 1] -> [- — -, oo[, x => G 2 ,*(x) = -. 

C C - 1 Cx - 1 

(39) 

With the minimal value for £ (Equation d37l i) we can slightly 
improve on the upper bounds for y+ and for vaj as 



,1/2 



C \ I/2 



-(cTl) ' ^-rtcTl) = Vph (cTl) 



(40) 

With the help of the information on the bounds for and £ 
we can refine the definitions given in Equations d2TT > and d22l > 
for F i ; * and Gi *, respectively, to their final versions as 

r C+l r n 1 / C \ 1/2 

^1* : . °°[-»] — . L 

u C- 1 L J V? \C + 1/ 



V2 

1 a 



+ 1 



V2\ Z 



1/2 



Gi,* : ] 



V2' 



C + 1 

c 



1/2 



y* => Gi,*(y*) 



] -»[ 



1 



. C+ 1 
C- 1' 



2y 2 * - l 



Let us now solve Equation < f30l > for x = 1/2R and find 
/ 1 £+1 



x = 



(41) 



(42) 



(43) 



2/? C £ - 1 

Equation d43l is a parametric representation of x = 1/2R in 
terms of f . As before we make this point explicit by introduc- 
ing the function F 2 ,+ defined as 

F 2 ,* : oo[-»R, f =» F 2 ,^)=i^-. (44) 



C- 1 



Cf - 1 



Since 



dF 2 



d£ C(£-l)< 



<0 



(45) 



Fi,*(0 is a decreasing function of In addition 

*2.*(fniin) = 1. Jim F 2 ,*(f) = p. (46) 

With this information on the function F 2j * we can refine its 
definition (Equation d44l) ). Combined with the definition for 
the function G 2 * (Equation d39b ) we obtain 

F 2 ,* : [|^|, 00 [ 1], £ => F 2 ,*(f) = 

(47) 

1 C + 1 Cx + 1 

G 2 ,* : ]-, 1] -» [— — , oo[, x => G 2 ,*« = 



C 



C- 1 



Cx - 1 
(48) 

F+2 is the inverse function of G 2 +, i.e., F 2> * = G^ and 
conversely G 2j * is the inverse function of F 2j *. 

4. ANALYTICAL SEISMOLOGY 
4.1. Summary 

Let us recapitulate the key results of the previous section. 
The three quantities that we assume to be known from obser- 
vations are the wavelength, A, the damping length, Ld, and 
the period, T. Alternatively, we can use the phase veloc- 
ity, v p h, instead of the period, T, if v p h is known instead of 
T. Analytical theory based on the TT and TB approxima- 
tions gives us two equations, namely Equations (0 and fl30"l > 
that express the wavelength A and the damping length Ld in 
terms of the density contrast, f , the normalized Alfven veloc- 
ity, = Vaj T/A = VA,i/v p h, and the inhomogeneity length 
scale normalized to the radius of the tube, x = 1/2R. These 
three quantities y+, and x are the seismic quantities in the 
sense that they are the quantities that we aim to determine with 
the use of observed values of the wavelength, A, the damping 
length, Ld, and either the period, T, or the phase velocity, v p h. 
Since we have only two equations that relate the three un- 
known quantities to the three observed quantities there are an 
infinite number of solutions. The seismic variables are con- 
strained to the following intervals 



X € I x - 



c + 1 



C-l 



C \ I/2 



J? \C + l) 



V2 



and are related to one another by 



y* = F i,*(f), ( = Gi,*Cy*), 
x = F % *(&, C = G 2 ,*(x) 



(49) 
(50) 
(51) 



(52) 
(53) 



The functions F l ir , G\ F 2 and G 2 ,+ are defined by Equa- 
tions (@J}, d42j, (@7J and |@8), respectively. 

In Equations d52l l and (l53l only two equations are indepen- 
dent since Gi,+ is the inverse function of Fi ^ and G 2>+ is the 
inverse function of F 2+ . Equations ( |52l i and ( T53l give us the 
infinitely many solutions of the seismic inversion in paramet- 
ric form. Each of the three unknowns can be used as parame- 
ter and the two remaining unknowns can be expressed in terms 
of that parameter. For example choose f as parameter. Let £ 
take on all values in 7^ and compute the corresponding val- 
ues of and x by the use of = F\^{Q and z = F 2j *(£). 
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Or choose y+ as parameter. Let y+ take on all values in 
and then compute the corresponding values of ( and x by the 
use of f = Gi,*(>*) and z = F 2 ,*(Gi,*(y*)). Finally, use x 
as parameter to define the solutions of the inversion problem. 
Let x take on all values in I x and then compute the corre- 
sponding values of y* and f by the use of f = G2,*(*) and 
y* = ^i,*(G 2 ,*0c)). 

Note that instead of the seismic variables y+, x, and the 
inversion can be performed for the variables v>A,i, l/R, and 
which have a more obvious physical meaning. The relation 
between both sets of seismic variables is 



A 
T 



VA,i = y+- = y*v ph , 



/ 

— = 2x, 
R 



(54) 



These alternative seismic variables are constrained to the in- 
tervals 



VA,i 



v P h / c y/ 2 

Vf' Vph lc7l] 



I 

— G 

R 



-l/R - 



Vph 

v! 

2 

c' 2 

C+ 1 



C- 1 



(55) 



(56) 



(57) 



4.2. Warning 

The present inversion scheme for propagating MHD waves, 
as its twin version for standing MHD waves, is simple to use. 
However, caution is required. First of all the analytic ex- 
pression used in the inversion scheme for the damping length 
or conversely for the damping time states that the damping 
length or damping time are inversely proportional to l/R. This 
relation is definitely very accurate when the non-uniform lay- 
ers are sufficiently thin. When these analytical expression 
are applied to fully non-uniform wave guides then the predic- 
tion would be that the waves undergo extremely fast damping. 
This result might be erroneous. Equation (l30t is derived for a 
thin non-uniform transitional layer and the kink MHD wave is 
essentially a surface Alfven wave. When the no n- uniform lay- 
ers become very thick the MHD waves are no longer surface 
wave s (see lArregui et al.l l2006t IVan Doorsselaere & Poedts 
120071) . Their damping deviates from that for surface Alfven 
waves and can be quite different from that predicted by Equa- 
tion (f30b . So it is wise to stay away from high values of l/R. 

Secondly we neglected non-uniformity along the wave 
guide. For standing MHD waves longitudinal non-uniformity 
in bo th density (see, e.g., lAndries et al.l 120051: lArregui et al.l 
2001 and magnetic fi eld (see, e. g.. IVerth & Erdelyil 120081: 



Ruderman et al.l 120081) affects the periods and the ratios of 



periods, but the damping by resonant absorpt ion is not af- 
fected by longitudinal stratification (see, e.g.. lAndries et alj 
120051: lArregui et alJl2005b iDymova & Rudermanll2006l) . Also 
longitudinal stratification does not cause additional damping 
or amplification. As the periods are concerned a longitudi- 
nally averaged density can be defined that produces the same 
period in a loop that is homogeneous in the longitudinal di- 
rection. For propagating waves the story is different. Radial 
stratification causes resonant absorption and damping but lon- 
gitudinal stratificatio n has an effect on the amplitude of the 
wave as shown by iSoler et all (1201 Id) . Longitudinal strat- 
ification causes the amplitude of the wave to increase with 
height and might partially or fully hide the damping due to 
resonant absorption. If the effect of longitudinal stratification 
is not removed from the observations, then Equation ( f30b uses 



and underestimate of the damping length due to resonant ab- 
sorption and this will affect the inversion result. 

Finally, we reca ll that the present invers ion scheme is based 
on linear theory. iRuderman et al.l (120101) showed that non- 
linearity can strongly increase the efficiency of damping due 
to resonant absorption. Although linear theory is accurate 
enough to describe t he small-amplitude w aves observed in 
the solar corona (e.g.. Tomczyk et al. 2007), the influence of 
nonlinearity may be important for kink waves propagating in 
other structures in the solar atmosphere. The conditions for 
which nonlinear effect s become important are discussed in 
IRuderman et all ( 120101) . 

5. EXAMPLE 

As an illustrative example of the technique described 
in the previous Sections we re-anal yze the CoMP obser- 
vations of running coronal w aves (|Tomczvk et al.l 120071 : 
ITomczvk & Mcintosh! 120091) . IVerth et alj (120101) showed 
that CoMP observations are consistent with an interpretation 
based on resonantly damped propagating kink waves. Here 
we perform the full seismological inversion using the CoMP 
data. 

We use VA,i, l/R, and £ as our seismic variables. To per- 
form the analytic inversion, we need observational values of 
the wavelength, A, the damping length, Ld, and the period, T. 
Alternatively, in the inversion of v>A,i we can use the phase 
velocity, v p h, instead of the period, T, because both quan- 
tities are related by the wavelength (Equation (TT4|)). Since 
th e observational value of v p h = 0.6 Mm s is provided 
by ITomczvk & Mcin tosh (2009J, we use v p h in the inversion 
scheme. 

Next we compute the parameter C = tiLq/A using observed 
quantities. In principle we need observational values of both 
the wavelength, A, and the damping length, Ld- However, we 
note that what we actually need is an observational estimation 
of the ratio L^/A. The theoretical expression of Lu/A is given 
in Equation ( f30b . Both A and Ld dep end on the wave p eriod, 
T, but their ratio is independent of T. IVerth et ail (2010) took 
advantage of this result and expressed the ratio Ld M in terms 
of the constant parameter £e as 



that allows us to rewrite C as 

A 



C 



(58) 



(59) 



The parameter ^e contains all the information about the prop- 
erties of the w aveguide. After p erforming a frequency fit to 
the CoMP data lVerth et alj ( 120101) found that the best estimate 
is £e = 2.69, wh i ch cor responds to C = 8.45. Note that the 
fit of IVerth etaD ((2010) was made fo r the wave path r epre- 
sented with a dashed line in Figure 1 of Tomc zyk & Mcintosh! 
(2009). The spatial resolution of CoMP observations was 
not enough to isolate individual coronal loops, i.e., indivi- 
tual waveguides, and the signal was spatially averaged in- 
si de the dotted regi on enclosing the dashed line (see details 
in IVerth et alfe OlO). Since the power was averaged inside the 
dotted region IVerth et alj (2010) obtained an averaged value 
of £e- For this reasion, the seismic variables inverted using 
this averaged value of ^e have to be interpreted as averaged 
values as well. T he estimated £e using the frequency fit in 
IVerth et ail d2010l) is consistent with previous estimations us- 
ing TRACE observations, which points out the validity of the 
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Figure 1. Result of the analytic inversion scheme applied to CoMP observa- 
tions of propagating coronal waves. The solid line is the ID solution in the 
3D space of parameters va,1-1/R-(- The dotted lines are projections of the 
full solution to the various planes. 

analysis. We stress again the ratio L^/A is enough for the in- 
version. However if information of both A and Lo is available 
in the observations, it can be used to directly compute C. 

Now we use Equations d55ll-(l57|i we compute the intervals 
of the seismic variables va;„ l/R, and namely 



v A;i e ]424, 567] km s" 1 , 
/ 



R 



e [0.24, 2], 



f e [1.27, oo[. 



(60) 

(61) 
(62) 



The variable that can be constrained in the narrower range is 
Va,i, whereas £ remains in practice unconstrained. The lower 
value of l/R is imposed by the observations while the upper 
value is imposed by the model. 

Now we perform the full inversion and use ( as a free vari- 
able. We compute the corresponding values of Vaj and l/R 
that are compatible with the observations. To do so we use 
the relations given in Equations (l52l l and (l53l l. namely 



VAi 



Vph /£ 



+ 1 



V2\ ( 



1/2 



/ 2 <r + 1 

-=2F 2 ,*(£) = - — . 



(63) 
(64) 



The corresponding solutions form a ID curve in the 3D space 
of parameters v^-l/R-^ (Figure[T|l. Any point on this curve 
is equally compatible with the observations. 

6. CONCLUSION 

In this paper we have presented an analytical seismologi- 
cal inversion scheme for propagating MHD waves in the solar 
atmosphere. This scheme uses the observational information 
on wavelengths and damping lengths in a consistent manner 
and is based on approximate asymptotic expressions for the 
theoretical values of both quantities. The scheme also needs 
observational values of the wave periods or, alternatively, of 
the phase velocities. The seismological inversion scheme for 



propagatin g waves shown here is the counterpart to that de- 
veloped by Goosse ns et alj d2008) for standing waves. 

Using the Alfven velocity, the inhomogeneity length scale, 
and the density contrast as seismic variables, we have shown 
that they can be constrained in intervals which depend on the 
observed values of wavelength, damping length, and phase 
velocity. The seismic variable that can be constrained the 
most is the Alfven velocity, while the density contrast is the 
less constrained variable. The whole collection of values of 
the seismic variables that are compatible with the observations 
describe a ID curve in the 3D space of variables. A priori, 
with no additional information on the uncertainties of the ob- 
served quantities, any point on this curve can equally explain 
the ob servations. As described in Arre gui & Asensio Ramos! 
( 201 ll) in the case of t he inversion scheme for standing waves 
(Goos sens et al.ll2008h . additional information can be consis- 
tently implemented in the Bayesian framework, which results 
in the full determination of the three unknowns, with correctly 
propagated uncertainties. As for standing waves, the use of 
Bayesian analysis with the present scheme for propagating 
waves is possible. 

The inversion scheme presented here is based on a sim- 
ple model for the magnetic wave guide. It is a challenge 
for future works to incorporate more realistic ingredients to 
the model in order to determine their impact on the seismo- 
logical inversion. In particular, the effec ts of flow (see, e.g., 
iTerradas et al.ll2010aHSoler et alj|201 lbl) . longitudinal strafi- 
ficat ion (see, e.g.,ISoleretaLll201ld) . and nonlinearity (see, 
e.g- lRuderman et al.l2010l) are worth being explored in forth- 
coming works. 
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